Bright and dark breathers in Fermi-Pasta-Ulam lattices 
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In this paper we study the existence and linear stabihty of bright and dark breathers in one- 
dimensional FPU lattices. On the one hand, we test the range of validity of a recent breathers 
existence proof [G. James, C. R. Acad. Sci. Paris, 332, Ser. 1, pp. 581 (2001)] using numerical 
computations. Approximate analytical expressions for small amplitude bright and dark breathers 
are found to fit very well exact numerical solutions even far from the top of the phonon band. 
On the other hand, we study numerically large amplitude breathers non predicted in the above 
cited reference. In particular, for a class of asymmetric FPU potentials we find an energy thresh- 
old for the existence of exact discrete breathers, which is a relatively unexplored phenomenon in 
one-dimensional lattices. Bright and dark breathers superposed on a uniformly stressed static con- 
figuration are also investigated. 
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I. INTRODUCTION AND MODEL 



Discrete breathers, also called intrinsic localized 
modes, are classical exact spatially localized time- 
periodic solutions which can be sustained by many non- 
linear lattices (see 0,01 for a review). In 1994 MacKay 
and Aubry Q rigorously proved their existence in Hamil- 
tonian lattices with anharmonic on-site potentials and 
weak coupling. Breathers arc obtained by continuation 
from the uncoupled case in which trivial breathers ex- 
ist (e.g. breathers with only one oscillator excited, the 
others being at rest). With the same technique, the ex- 
istence of breathers was alsoproved for diatomic Fermi- 
Pasta-Ulam (FPU) chains In this model, two dif- 
ferent masses alternate on the chain and are nonlinearly 
coupled to their nearest-neighbours via an interaction po- 
tential V. This result is valid for a large mass ratio, 
since light masses (next-nearest neighbours) are weakly 
coupled due to the presence of a heavy mass in between. 

Unfortunately this method is not applicable to homo- 
geneous FPU lattices, which do not possess an uncoupled 
limit in which trivial breathers exist. For some years the 
only existence result concerned the particular potentials 



Vix) 



with m > 2 5] . FPU lattices seemed elusive 



to a rigorous mathematical treatment in spite of numer- 
ous approximate ^ J^, iSj and numerical ^ 0| studies, 
which have indicated the existence of discrete breathers 
in these systems. 

Nevertheless, recent papers have presented rigorous ex- 
istence proofs for discrete breathers in infinite FPU lat- 
tices. On the one hand, Aubry et al have proved 
the existence of breathers with frequencies above the 
phonon spectrum, when U is a strictly convex polyno- 
mial of degree 4. These results are obtained via a vari- 
ational method and apply in fact to higher dimensional 
generalizations of FPU lattices. Without additional as- 
sumptions on V, these results give only a partial infor- 



mation on breathers amplitudes. Under the additional 
condition that V is even, Aubry et al prove the exis- 
tence of breathers of arbitrarily small amplitudes in one- 
dimensional FPU lattices. 

On the other hand, the existence (resp. non ex- 
istence) of small amplitude breathers with frequencies 
slightly above the phonon band has been established 
when V satisfies (resp. violates) a local hardening condi- 
tion The proof is based on a centre manifold tech- 
nique, which shows more generally that all small am- 
plitude time-periodic solutions of the FPU system (in- 
cluding breathers) are determined by a two-dimensional 
map, provided their frequency lies near the phonon band 
edge. The discrete centre manifold method has been put 
in a general framework and applied to other systems, 
namely diatomic FPU chains (far from the uncoupled 
regime) and spin lattices [T5l |. 

The aim of this work is to test numerically the range of 
validity of the centre manifold method jl^ and to explore 
new phenomena, far from the small amplitude regime. 

The one-dimensional FPU system is given by the equa- 
tions: 



in = V'{Xn+l - Xn) - V (Xn 



-i), nez, (1) 



where a;„ represents mass displacements from their equi- 
librium positions and U is a smooth interaction potential 
satisfying V'{0) = 0, U"(0) = 1. In ref. [H, the existence 
of small amplitude breathers (SAB) with frequencies uJb 
slightly above the phonon band {uJb > 2) is obtained for 
B > 0, where 



S = lu(4)(0)-(U(3)(o))2 



(2) 



and their non-existence is proved for B < (see |Tg, 
17] for related results on the modulational instability of 
nonlinear normal modes and the tangent bifurcation of 
standing waves respectively). 
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The parameter B can be interpreted as a hardening 
coefficient, since breathers with amplitude A « have 
frequency Wf, ~ 2 + ^ A^. Note that B is slightly different 
from the classical hardening coefficient of an anharmonic 
potential {V is hard if f V^") (0) - (V"^^) (0))^ > and soft 
for |l/(4)(0)-(l/(3)(0))2 <0). 

In this paper we consider anharmonic potentials 



K- 



(3) 



(at least one of the coefficients K^, K4 is nonzero). One 
has B = 3K4 — AK^ and thus the parameter B is positive 
if jivTsI < \/3K4/2. In the case K3 = we have an even 
potential and the sign of B coincides with the sign of the 
quartic coefficient K4. 

In section II we compute numerically SAB and find 
that they are very well fitted by leading order analyti- 
cal expressions deduced from ref. j^^l- section III we 
consider a class of potentials for which numerical compu- 
tations yield large amplitude breathers (LAB) with fre- 
quencies near the top of the phonon band. In this case 
there is an energy threshold for the existence of exact 
discrete breathers. The existence of an energy threshold 
for breather creation in higher-dimensional lattices is well 
known 18, 19], but only a few one-dimensional examples 
have been given where such a threshold exists plll20ll2ll| . 
Section III explores in more details the example of ref. 
[2lj |. Section IV deals with dark breathers, i.e. spatially 
modulated standing waves whose amplitude is constant 
at infinity and vanishes at the center of the chain. Lead- 
ing order analytical expressions for small amplitude dark 
breathers (with frequencies inside the phonon band) fit 
the numerically computed dark breathers very well in the 
case of even potentials. In section V we analyze bright 
and dark breathers superposed on uniformly stretched 
or compressed static states. In each section the linear 
stability of the above mentioned solutions is numerically 
investigated. 



The main results of [13| can be formulated as follows. For 
ujh fa 2 (maximal phonon frequency), all small amplitude 
solutions Un (even in t) have the form 



y-nit) = bnCOs{uJbt) + ip{t, b^^i , bn, UJb) , 



(7) 



where (p is a smooth time-periodic function with 
frequency LUb and Sup | (^(i, 6„_i, 5„, cj^) | = 

0((|6„-i| + K\){K-i\ + K\ + K - 2|2)). This re- 
suit can be seen as an exact version of the rotating-wave 
approximation method |^, in the small amplitude limit. 
Note that ip can be computed at an arbitrary order (see 
[T^j l p. 51). Moreover, &„ satisfies the second order non- 
linear recurrence relation 



26„ 



= —fi bn + B + h.o.t 



(8) 



where /i 



-4^1 and the coefficient B is given in jSJ. 



Higher order terms are 0( |6„| ((|6„_i| -I- |5„|) -I- |^|) ) 
(note that the right side of lO can be computed at an 
arbitrary order). Consequently, the problem of find- 
ing small amplitude time-periodic solutions of © re- 
duces to the problem (|SJ), which can be viewed as a two- 
dimensional map {bmbn^i) — > (6„+i,6„). 

For B > Q fixed and /i > sufficiently small, the recur- 
rence relation has homoclinic solutions to satisfying 
lim 6„ = (a proof has been given in [13], using 

n — ^±cxD 

the invariance n — > —noi^). More precisely, there exist 
homoclinic solutions to denoted as ±&^, ±6,^, having 
different symmetries 6l_„_i = —b\, btn = 

Beyond this particular result, one can expect the ex- 
istence of infinitely many homoclinic solutions to (not 
necessarily symmetric). Indeed, for B > the intersec- 
tions of the stable and unstable manifolds of 6„ = are 
generically transverse, which yields the existence of "ho- 
moclinic tangles" p^ . 

A formal way to understand why homoclinic orbits ex- 
ist is the following. Setting ^„ = (— 1)"6„ yields the 
recurrence relation 



II. SMALL AMPLITUDE BREATHERS 

A. Summary of the theory 

We consider time-periodic solutions of (QJ having a 
given frequency Wf,. We introduce the interaction force 



where 



(4) 



(5) 



The time-average of y„ is independent of n (integrate 
(0), and is fixed to in reference T^, as it is the case 
for all spatially localized solutions. Then, problem 
leads to 



{y') ^iUn) = Vn+l - 2y„ -I- Vn-l, 71 G Z. (6) 



Cn+1 + ^ 2^„ — fi^n ^ B£^„ 



h.o.t. 



(9) 



Since < fi <^ 1, equation ^ can be approximated by 
an integrable differential equation 



with 



u(n07). 



(10) 



(11) 



Equation (|10|l has the homoclinic solutions v{x) — 
±V2/ cosh(a; + c). 

Using the above analysis, it has been proved that for 
B > and Wf, — 2 > sufficiently small, there exist SAB 
solutions of lO given b y eq uation {Tj) and homoclinic so- 
lutions of ||SJ) (see |l3 , ) • SAB have the form 



Un = (-l)"Cn COs{uJbt) + h.o.t. 



(12) 
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where £,„ = satisfies the recurrence relation @ 

and l^nl ^ as n — > ±00. The above homochnic so- 
lutions 6^, 6^ correspond via equation to SAB y^, 
having different symmetries j/^„_i(t) = yl^{t + Ti, / 2) , 
= y^, where Tj, = 27t/uji, is the breather period (see 
Pl l p.58,59). Note that y„(i + Th/2) ^ -ynit) in general 
due to higher order terms in (|12|l (however the egality 
holds if V is even). Homoclinic solutions — fojj of © sim- 
ply correspond to yl^{t + / 2) . 

The exact solutions y,\, can be approximated at 



leading order using l(Tn|l - l(TT)l with v{x) 
This yields 



\/2/ cosh: 



'2/i coswfji 



(13) 



S cosh(ny//i) ' 
and by symmetry we construct the second approximation 



yn{t) (- 



COS ujbt 



B cosh((|n + l/2|-l/2)^)' 



(14) 



Under these approximations one has yl^{t + Tjy/2) « 
—yliit) and ?/o(i) ~ yo(^)- ^^^t, the exact solutions 
satisfy yl{t + n/2) = -yl{t) + 0(|^|) and %i(0 = 
2/0^(0 +0(|/ip/2) as/i^O+. 

Notice that approximation p3|) can also be derived us- 
ing multiscale expansions [l^ . Moreover, these expres- 
sions also approximate w„ since Wn = ?/« + ^(y,^ ). 

These calculations show that the maximum amplitude 
Aw y^of the SAB is 0(|tJ6-2|i/2) asw6-^2+. More 

precisely, breathers with amplitude A sa have frequency 
u;t,-2+^A\ 

Although SAB decay exponentially as \n\ +00, their 
width is 0(|cjf, — 2|^^/^) and diverges as Wh 2+. As a 
consequence, if i? > breathers exist for any small value 
of energy in the FPU system . 



B. Numerical continuation 

We have computed these solutions numerically in order 
to check the range of validity of approximations (|13|l . 
(fn|l . We have fixed ^4 = 1 in eq. Q and thus the 
parameter B is positive if \K^\ < — \/3/2 ~ 0.86. 

We have performed our computations using a numeri- 
cal scheme based on the concept of anti-continuous limit 
and path continuation with the Newton method |23j . 
With this technique it is more convenient to use the dif- 
ference displacements variables u„. Indeed, with these 
new variables the dynamical equations become 

iln + 2V'{Un) - C [V'{u„ + i) + V'iUn^l)] = , n G Z , 

(15) 

where C = 1 for our system, but it makes possible to 
consider the system ()15|l as a Klein-Gordon system, with 
an on-site potential V and an anticontinuous limit at 
C = 0. Note that u„ is one-to-one related to the forces 
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FIG. 1: Comparison between an exact site-centred mode (cir- 
cles) and the approximation (1131 (dashed line) for uib ~ 2.01 
{jj, ~ 0.04). The cubic coefficient of the potential is K3 = —0.3 
{B — 2.64), but similar results are obtained for any other fixed 
value of K3 with {K^l < {B > 0) and jj. > small enough. 



yn at small amplitudes, because V' is locally invertible 
since ^"(0) ^ 0. 

We use periodic boundary conditions itn+2p(0 — ^n{t)- 
The periodicity is considered for an even number of 
sites in order that the maximum frequency of the lin- 
ear phonons is exactly 2 as in the infinite lattice (this 
frequency corresponds to the phonon with wave number 

We have obtained breathers with symmetries Un(t) — 
w_„(t) (site-centred mode) and u„(i) = M_„_i(t -I- Tb/2) 
(bond-centred mode). Note that the coordinates trans- 
formation (O produces an exchange of the symmetry 
properties between both modes, i.e. in the difference dis- 
placements variables the site-centred mode is the Page 
mode and the bond-centred mode is the Sievers- 
Takeno mode 0. 

FigQ] shows the excellent agreement between the ap- 
proximation (|13(l (dashed line) and an exact site-centred 
mode obtained numerically (circles) for uji, = 2.01 (/z « 
0.04). 

We have continued the SAB obtained numerically for 
B > as cjfc goes away from the phonon band. In 
Fig|21 we have plotted the maximum amplitude of the 
forces sup|j/„| (dash lines) versus /i^/^ for two different 
values of the cubic coefficient, K3 = —0.6 {B ~ 1.56) 
and K3 = —0.3 {B = 2.64) from top to bottom, re- 
spectively. The maximum amplitude of the forces is an 
approximately linear function of /.t^/^ up to /i « 0.6. The 
continuous lines represent the maxima given by approx- 
imation (|13|l . This approximation is not satisfactory for 
very small values of B since higher order terms in © 
cannot be neglected when B ^ 0. More precisely, the 
smaller B is, the smaller fi must be chosen in (|13|) in 
order to have a good approximation. 

The maximum amplitude of the relative displacements 
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FIG. 2: Force (dashed line) and relative displacements 
(squares) maxima versus The continuous line corre- 

sponds to the maxima given by approximation 11311 . Pa- 
rameters from top to bottom: — —0.6 {B = 1.56) and 
Ks = -0.3 {B = 2.64). 



FIG. 4: Maximum amplitudes, sup\un\, of the Sievers-Takeno 
mode versus frequency for an even potential. The continuous 
line is the maxima given by eq. CZJ, while the dashed line 
corresponds to the Sievers-Takeno approximation. Squares 
correspond to the numerically computed solution. 



T 1 1 r 



1 

0.5 

Un(0) cteoee-eoeooo-H- /' 

b 

-0.5 
-1 



1 9 



J I I L. 



_i_ 



-15 -10 -5 5 10 15 
site 



FIG. 3: Comparison between an exact bond-centered mode 
(circles) and the analytical approximation (1141 (dashed line) 
for jj, = 1.76. The value of fj, corresponds to a breather fre- 
quency ijJi, — 2.4 relatively far from the phonon band. The 
cubic coefficient of the potential V is Ks = —0.3 {B = 2.64). 



sup|m„| is also an approximately linear function of /i^/^ 
(see squares in fig|21l. This was expected for small fi 
since Un = yn + O(y^), but it occurs far from the phonon 
band, at least until values of /i w 2.25 . In fact, we have 
checked that, surprisingly, expressions (|13|l and H14|l fit 
very well the profile of the relative displacements u„ even 
far from the phonon band as Fig. Oshows. Note that the 
vibration amplitudes of the breather in Fig. |31are quite 
large. 

In Fig 2] we have plotted again the maximum ampli- 
tude of the relative displacements (squares), but now for 
an even potential [K^ — 0) and larger frequencies, in 
order to compare approximation (|14|l with an approxi- 



mation derived by Sievers and Takeno in the large am- 
plitude regime ^ . The Sievers-Takeno approximation is 
obtained using the rotating wave approximation (only the 
first Fourier component of solutions is taken into account) 
and Green's function techniques. For the Sievers-Takeno 
mode and = 1 {B = 3) it predicts 



SUp\Un\ 



(16) 



This amplitude-frequency relation corresponds to the 
dashed line in Fig 01 On the other hand, from eq. H14|) 
it follows (continuous line) 



SUp\Un\ 



(17) 



Clearly we see that approximation (|16|) fits better very lo- 
calized solutions with large amplitudes, whereas approx- 
imation (|17|l works better up to moderate amplitudes 
{ujb « 2.4 in this case). In this sense both approxima- 
tions can be considered complementary. 

In addition we have performed a Floquet analysis in 
order to study the linear stability of the SAB |23|. In 
the symmetric potential case the result is well known : 
the Page mode is stable while the Sievers-Takeno mode 
has a harmonic instability (a pair of real eigenvalues 
CT, (T~^ close to 1) that increases with the breather fre- 
quency. When a cubic term in the potential V is intro- 
duced {K-i 7^ 0), the situation is more complex. The 
Sievers-Takeno mode shows again a harmonic instability 
(see FigEt) but it also appears oscillatory instabilities. 
As FigO) shows for the Page mode, these oscillatory in- 
stabilities increase with \K^\. Although most of them are 
size dependent effects due to the discretization of contin- 
uous spectrum psj , for the infinite system some of them 
might remain out of the unit circle and in this case SAB 
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FIG. 5: a) Harmonic instability of the Sievers-Takeno mode 
for Kz — —0.3; b) Oscillatory instabilities of the Page mode 
versus K3 for a lattice with 52 sites (full circles) and a lattice 
with 100 sites (pluses). The frequency is Wb — 2.1. 



would be unstable. We shall not examine this problem 
here and leave it for future works. 



III. LARGE AMPLITUDE BREATHERS WITH 
AN ENERGY THRESHOLD 

In this section we study breathers in the parame- 
ter region B < 0. Using the recurrence relation JHJ, 
one can show that SAB with frequencies slightly above 
the phonon band do not exist for i? < (see [T^ . 
[Tsf). However, large amplitude breathers (LAB) do ex- 
ist in this parameter region |^. More precisely, for V 
strictly convex {K4 > 0, IK^] < y/SKl) and B < 
{^/3Kl/2 < IK^I), breathers whose amplitudes do not 
tend to zero as Wb — > 2+ have been numerically obtained 
a consequence there is an energy thresh- 
old for breathers creation in these FPU systems. The 
criterion ^/3K4/2 < {K^l < \/3Kl for the existence of 
an energy threshold originates from reference [2]| . and 
follows by combining analytical results in (non exis- 
tence of small amplitude breathers) and [lll | (existence 
of large amplitude breathers). 



Un(0) 




FIG. 6: Comparison between a small amplitude breather (full 
circles) for B > and a large amplitude breather (blank 
squares) for B < with the same frequency wt — 2.01. 
Parameters: SAB: K3 = -0.3, K4 = 1 {B ^ 2.64); LAB: 
Kz = -1, K4 = 1{B = -1). 



The existence of an energy threshold for breather cre- 
ation has been formally analyzed by Flach, Kladko and 
MacKay for general Hamiltonian systems, in a dif- 
ferent case when small amplitude breathers bifurcate at 
an edge of the phonon spectrum. An energy thresh- 
old exists above a critical lattice dimension d (typically 
d = 2) which depends on the nonlinear terms of the 
Hamiltonian. This analysis has been extended to a class 
of non-analytic Hamiltonian systems in fj^l- The ex- 
istence of excitation thresholds has been rigorously an- 
alyzed by Weinstein il^J for the d-dimcnsional discrete 
non linear Schrodinger equation, using variational meth- 
ods (see also for formal variational arguments in the 
one-dimensional case). 

In Fig. |n|we compare a SAB and a LAB with the same 
frequency Wfc = 2.01. As it is shown in Fig [7^ the energy 
of a LAB (dash line) does not go to zero when approach- 
ing the phonon band edge. It lies always above a certain 
positive lower bound. In the case of SAB (continuous 
lines) the breather energy can be lowered to arbitrarily 
small values. 

The dependence of the energy threshold on the asym- 
metry parameter is also shown in Fig. [TJa, where we 
have fixed K4 = 1. One observes that the energy thresh- 
old goes to as B ^ 0" (for K3 -\/3/2 « -0,86). 
Moreover, LAB families also exist for nonconvex poten- 
tials {K3 < —VS « —1,73), and their energy thresh- 
old increases with \K3\. This shows that the criterion 
V 3K4 1 2 < jA'al < \fTK4 for the existence of an energy 
threshold is not optimal. We conjecturate that an en- 
ergy threshold appears under the more general condition 
K4 > 0, V3A'4/2 < lATsI {B < 0) in the one-dimensional 
FPU system. 

Figure IHl shows the energy of a LAB (upper line) and 
its amplitude (lower line) versus frequency, close to the 
phonon band edge (the upper line is a close-up of figure 
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FIG. 7: (a) Breather energy versus frequency close to the 
phonon band edge (we fix K4 = 1). The dashed hue shows 
the existence of a threshold in energy for a LAB with K3 — ~1 
{B < 0). The continuous lines {K3 = —0.6 and Ks — from 
top to bottom) show that for B > SAB of arbitrary low 
energy exist, (b) Dependence of the LAB energy threshold 
on the asymmetry parameter K3. 



near Ub — 2). The reason for the non monotonous 
behaviour of energy is that there are two competing ten- 
dencies: on the one hand amplitudes are decreasing with 
frequency and going to a nonzero lower bound, but on 
the other hand the width of the breather is increasing 
and more and more particles begin to oscillate as ujh de- 
creases. 

We have found LAB with the same symmetries as SAB 
(Page and Sievers-Takeno modes). LAB have an expo- 
nential decay for cuf, > 2. As Fig. |51 shows, we have been 
able to fit the LAB profiles to the following expression in 
the case of the Page-mode: 



Un{t) « Q;(wfc) o-(a;h)l"l cos (ojbt) 



(18) 



where 



Wf, 2 



a(c.fc) = l-:^ + y(c.,^-4)i/2e(-l,0) , (19) 

and a{ujb) is a fitting parameter. A similar expression 
corresponds to the Sievers-Takeno mode. 
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FIG. 8: LAB energy (upper line) and amplitude (lower line) 
vs frequency close to the phonon band edge (we fix K3 = — 1, 
Ki = 1). 



We have computed the decay rate |tT(a;(,)| as follows. 
The FPU system ^ can be formulated as a mapping 



Vn 

Vn+l 



= F 



Vn-l 



where the map 



F 



Vn-l 
Vn 



{V')~Hyn) + 2yn-y„- 



(20) 



(21) 



acts in a space of smooth time-periodic functions with 
frequency > 2. The linearized operator 



i?F(o)r""M = 

V Vn J 



Vn 

2)y„ 



Vn-l 



(22) 



has a purely hyperbolic spectrum and cr(wh), (T(ajfc)^^ are 
the closest eigenvalues to —1 (with a (2) — —1). 

We have observed that LAB remain localized for oji, — 
2 (see pluses in figure O and thus approximation l|18|) 
does not work in this limit. 

The stability properties of LAB are the same as SAB 
with non-even potentials, i.e. both modes show oscilla- 
tory instabilities and also a harmonic instability in the 
case of the Sievers-Takeno mode. As we mention in sec- 
tion II, most of the oscillatory instabilities are size de- 
pendent effects due to the discretization of continuous 
spectrum "2^. However, a band analysis 0| suggests 
that, at least for some parameter values (close to = 3, 
A'3 = — 1, A'4 = 1), some eigenvalues remain out of the 
unit circle for the infinite system and LAB are unstable. 
We have studied the evolution of Floquet eigenvalues of 
the Page mode for uj^ ~ 3, if 3 = — 1 and K4 — 1 when 
the system size becomes larger. In figure [101 one can see 
that one pair of eigenvalues seems to converge to values 
out of the unit circle for large system sizes, as stated in 

m 
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ln|u,(0)| 




general different from the zone-boundary modes a;„ (t) = 
(— (also called binary oscillations). Indeed, in 
this latter case the time-average of relative displacements 
Un is 0, whereas interaction forces have time-average in 
our case. However, both types of solutions coincide for 
even potentials V. 

To obtain an approximate expression of the hetero- 
clinic solutions, we can rewrite the recurrence relation 
in the following form 



Pn+l + Pn-1 — 2p„ 



Bpl + h.o.t. 



(23) 



using the variable change ^„ = y/—p p„. Equation H23|) 
can be approximated by the differential equation: 



FIG. 9: Logarithms of the Page mode profiles for Kg = —1, 
K4 = 1 and Wb — 3 (circles), uit = 2.1 (squares), wt = 2.05 
(triangles) and Wb — 2 (pluses) . The dashed line is the 
logarithm of eg. (1181 after fitting the parameter a. Note that 
eg. 1181 does not work in the limit cub ~ 2. 



v" = -v- Bv-^ 



(24) 



with pn = v{ny/JI). Consequently, the system has a 
family of small amplitude solutions which can be approx- 
imated by 
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FIG. 10: Floquet eigenvalues of the Page mode for uob — 
3, ifs = —1 and K4, — 1, depending on system size. The 
eigenvalue of largest modulus and its inverse seem to converge 
to values out of the unit circle for large system sizes. 



IV. DARK BREATHERS 

Using the recurrence relation JSJ, it has been shown 
that for S < and /i < small enough (w;, lies in the 
phonon band) there exist small amplitude heteroclinic 
solutions of connecting two nonlinear normal modes 
at infinity, i.e. liin yn{t)~y^{t) = (see [T^ . 



n — ^±00 



We now briefly describe the nonlinear normal modes 
y^. These solutions are spatially periodic iyt+2(t) = 
yn{t)), have time-average (the mean interaction force 
is 0) and neighbours oscillate out of phase (yt+ii^) — 
y^{t + Tb/2)). Moreover, solutions y^ are equal up to 
a half period phase shift {yn{t) = 2/n + ^b/2))- Note 
that the corresponding displacement patterns Xn are in 



yn{t) 




cos{uJbt)- (25) 



This type of solutions are called dark breathers [23| or 
antisymmetric kinks Typically the term breathers 

refers to bright breathers, that is a few excited units in 
a lattice. Here we have the opposite kind of localization: 
most oscillators are excited except one or a few of them. 

Expression (|25|l approximates exact site-centred solu- 
tions of ((nj satisfying y-n{t) = yn{t + Tb/2). Under our 
approximation one has yn(t+Tb/2) ^ ^yn{t), and in fact 
the exact solutions satisfy yn{t + Tb/2) = —yn{t) + 0{\p,\) 
(uniformly in n € Z) as /i ^ 0^. 

Equation (jSJ has also exact bond-centred dark 
breather solutions satisfying 2/„(t) = 2/_„+i(i). We shall 
approximate them by 



yn{t) ~ (- 




(-n+ 1/2)7^ 



cos{ujbt)- 



expressions 
Vn + Oiyl). 



V2 

(26) 

(j2SJi l|26|l also approximate Un 



Note that 
since u„ 

For computing dark breathers numerically we consider 
a lattice with an odd number of sites (51 lattice sites). 
We first choose an even potential {K^ = 0) with K4 < 
(so that B < 0). We solve equation (|15|l in Fourier space 
by the Newton method, fixing the time-average of uq to 
and using (|25() as an initial condition. Dark breathers are 
computed for non-even potentials by continuation with 
respect to K3, starting from = (during the contin- 
uation the time-average of y„ remains very small, due to 
the constraint on uq). 

Fig.^Jshows the excellent agreement between an exact 
site-centred solution obtained numerically (circles) and 
the approximation (|25|) (dashed line) , for an even poten- 
tial with K4 = —1 {B = —3) and a frequency ujb — 1-9. 
For non even potentials (full circles) the agreement is not 
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FIG. 11: Small amplitude site-centred dark breathers for a 
frequency u>b = 1-9. The circles represent an exact solution 
computed numerically with an even potential (K/t = —1), 
while the full circles were obtained adding to the previous 
potential a cubic term with a coefficient Kz = 0.6. These 
solutions have the symmetry u_„(t) = Un(t ^r Tb/2). The 
dashed lines are the tanh function of approximation 125^ for 
both cases. 



FIG. 13: Maximum amplitudes, sup|u„|, of the site-centred 
dark breather versus —^^^"^ for Kz ~ (circles) and Kz = 0.6 
(squares). In both cases if 4 = —1. The dashed line is the 
maxima given by eq. 12511 . In the non-even case approxima- 
tion H25|l does not provide a good fitting far from the phonon 
band edge. 
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FIG. 12: Small amplitude bond-centred dark breather for 
a even potential with K4, = — 1 and a frequency uJb = 1.99. 
The circles represent an exact solution computed numerically, 
having the symmetry u„(t) — it_„+i(t). The dashed line 
corresponds to the analytical approximation H26|l . 



SO good. The reason is clearly that the approximation of 
Un by Un with equation l|25(l has the symmetry property 
Un{t + Th/2) ~ —Un(t) which is not true any more when 

Similar results are obtained for bond-centred solutions 
(see Fig.ini. 

For a fixed potential we can continue numerically these 
solutions decreasing uJb- As Fig. 1131 shows, then the os- 
cillation amplitudes increase and for even potentials ex- 
pression (|25|l continues matching very well the relative 
displacements u„ (circles), even far from the top of the 
phonon band. For non even potentials (squares) the 
agreement is also excellent close enough to the edge of 



the phonon band. For lower frequencies we observe dis- 
crepancies becoming larger. 

We have also performed a Floquet analysis to investi- 
gate the linear stability of these family of solutions when 
we decrease the frequency. 

On the one hand, the site-centered mode is stable close 
to the edge of the phonon band for even potentials as 
figure ink shows. However, a cascade of harmonic insta- 
bilities appears for low enough frequencies (larger am- 
plitudes). These instabilities are visible in Fig. 114b at 
u!b ~ 1.43. Before coming out the harmonic instabilities 
one observes an instability bubble at tob ~ 1.53. Nev- 
ertheless, its magnitude decreases as the system size in- 
creases and thus it should not be relevant for the infinite 
system. 

On the other hand, the bond-centered mode shows a 
harmonic instability (see figure IT^) that increases when 
the frequency decreases. One also observes another cas- 
cade of harmonic instabilities at LOb ~ 1.43. 

Non-even potentials induce oscillatory instabilities for 
low enough frequencies (larger amplitudes), both for the 
site-centered and the bond-centered modes. Fig. I^b 
and fig. 115b show that this occurs for frequencies lower 
than 1.85, for K3 = —0.6 and K4 = —1. Numerically 
we haven't found these oscillatory instabilities near the 
phonon band edge, but we do not exclude that they might 
exist and be harder to detect. 

The harmonic instability before mentioned for the 
bond-centered mode is also present for non even poten- 
tials, but is not visible in fig. 115b due to oscillatory in- 
stabilities. Note that frequency scales in figures a) and 
b) are different. For this reason in figures b) we do not 
see either the cascade of harmonic instabilities around 
UJb ~ 1.43. We have checked that they are also present 
with non-even potentials. 
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FIG. 14: Evolution of the Floquet eigenvalue moduli versus 
frequency for the site-centred dark breather, a) Even poten- 
tial {Ks — 0); b) Non-even potential {K3 = —0.6). In both 
cases if 4 = — 1. 
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FIG. 15: Evolution of the Floquet eigenvalue moduli versus 
frequency for the bond-centred dark breather, a) Even poten- 
tial {Ks ~ 0); b) Non-even potential {K3 — —0.6). In both 
cases K4 = — 1. 



V. BRIGHT AND DARK BREATHERS WITH 
AN UNIFORM STRESS 

Breathers considered in section III for B < can be 
numerically continued at fixed frequency up to values 



B > 0, bringing about a new type of solutions which 
satisfy lim„^±cxD Un = c, where c is a nonzero constant. 
As Fig. [TCI shows, these new type of solutions can coexist 
with SAB for the same values of the parameters and can 
be seen as time-periodic oscillations around static solu- 
tions of ^ given by The lattice is uniformly 
stretched at infinity for c > 0, and uniformly compressed 
for c < 0. Such breather solutions with a constant static 
strain c < have been also observed by Sandusky and 
Page in the case B <0 (see p.876). 

This section provides analytical and numerical re- 
sults concerning breather solutions with a constant static 
strain. The case of dark breathers with an uniform stress 
is also examined. 

In the sequel we assume V"{c) > 0. We make the 
change of variable — cn + Xn {t^/W{cj) where 

i„ <C 1 and introduce the modified potential 



V{u) = {V{c + u) - V'{c)u)/V"{c). 
Equation ^ leads to 



(27) 



V'{Xn+l-in)'V'{x n ^n — l 



^eZ, (28) 



with V'{0) ~ and ^"(0) = 1. One can see that equation 
(|28|1 is exactly equation where V has been changed to 
V. Consequently, the results of the above sections (both 
analytical and numerical) readily apply to equation l|28f) . 

Returning to the original variables, the following an- 
alytical results for the FPU system |^ follow from the 
analysis of section II-A. Let us fix c e M and define 



Bic) = ^V"{c)V^''Hc)-iV^'Hc)r. 



(29) 



If B{c) > and - 2y/V"{c) > is sufficiently small, 
the FPU system has a family of exact solutions which 
can be approximated by 



u„{t) 



2)1 



cos(ujbt) 



B{c) cosh{ny/fi/V"ic)) ' 



(30) 



where fi — uj^ — AV" {c), < ^ 1. These solutions are 
time-periodic (with frequency tUb) and their oscillating 
part is spatially localized. As n ^ ±oo they converge 
towards a uniformly stressed static state (stretched or 
compressed depending on c). Note that their lower fre- 
quency lies inside the phonon band if V"{c) < 1, and 
above for V"{c) > 1. The exact solutions approximated 
by (|30|) are site-centred, i.e. they satisfy Un{t) = U-n{t). 
There exists also exact bond-centred solutions satisfying 
Un{t) = u-n-i{t + Tr/ojb). In the sequel we shall denote 
these solutions as breathers with an uniform stress. 

As we previously mentioned, the numerical study of 
these solutions is already contained in section II-B, which 
applies to the renormalized displacements i„ satisfying 
(|28|l . For the polynomial potential breathers with an 
uniform stress have the same stability properties as the 
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FIG. 16: Comparison between a SAB (full circles) and a 
breather with an uniform stress (continuous line) for the same 
values of the parameters: Kz — —0.6, K4 = 1, wt = 2.01. 

SAB corresponding to the modified potential H27() . In 
particular, oscillatory instabilities are observed for K3 — 
{V is even), since the modified potential V is non-even 
in this case. On the contrary, if K3 ^ 0, K4 > if|/6 
then the modified potential V is even for c = —-^j^ (and 
the conditions V"{c) > 0, B{c) > are realized). In this 
case oscillatory instabilities are not observed for the SAB 
with specific static strain c. 

If B{c) < and K4 > 0, we numerically obtain 
large amplitude solutions of the same type, whose os- 
cillating parts do not vanish as Wb 2y/V"{c) (with 
Wb > 2y^V"{c) ). This result follows directly from the 
numerical study of section 111, as well as the stability 
properties of these solutions. These LAB with an uni- 
form stress have been previously observed by Sandusky 
and Page (see p. 874-876), for the interaction poten- 
tial © with B < 0. In their examples V"{c) > 1 (since 
c < 0, K3 < 0, K4 > 0) and thus the lower frequen- 
cies of these LAB families lie strictly above the phonon 
band. This phenomenon can be observed in fig. 8 of ref- 
erence p^, which provides frequency vs amplitude plots 
of such solutions (note that c may not be constant on 
each curve). 

If B{c) < and uib - 2y/V"{c) < is sufficiently small, 
the FPU system has a family of exact solutions which 
can be approximated by 

.„(.) . c + (-D-y^tanh (^^) 

(31) 

where ^ = - 4F"(c) < 0, |^| < 1. These solu- 
tions are time-periodic (with frequency Ub) and consists 
in a dark breather superposed with a uniformly stressed 
static state. Note that their higher frequency lies above 
the phonon band if V"{c) > 1. These solutions are site- 



centred {u-n{t) = Un{t + IT /iUb)) , but their bond-centred 
analogues also exist. They have the same stability prop- 
erties as dark breathers in section IV, for the modified 
potential (|77|l . 

VI. CONCLUSIONS 

In this work wc have analyzed both numerically and 
analytically a rich variety of nonlinear solutions in FPU 
lattices. There are four main results. 

On the one hand, we have numerically explored the 
range of validity of a recent breathers (bright or dark) ex- 
istence proof [121 for one-dimensional FPU lattices. Our 
computations have been carried out with polynomial in- 
teraction potentials of degree 4. For i? > 0, we have 
found that the maximum amplitude of bright breathers 
in the relative displacement variables is an approximately 
linear function of /i^^^ even far from the phonon band 
and, in fact, the approximate expression (|14|l describes 
very well the profile of the breathers when their ampli- 
tudes are large. For B < and K4 < 0, we find in the 
same way that dark breathers can be very well described 
by approximate analytical expressions. 

On the other hand, we have found numerically for 
B < large amplitude bright breathers having frequen- 
cies arbitrary close to the top of the phonon band, which 
are out of range of the local analysis ■ This breather 
family exhibits an energy threshold, which is a rarely ob- 
served phenomenon in one-dimensional lattices. 

Moreover, with respect to stability properties we have 
checked that non-even potentials induce oscillatory in- 
stabilities in both bright and dark cases. In the case 
of even potentials, the Page mode is linearly stable while 
the Sievers-Takeno mode is unstable. We have also found 
that bond-centred dark breathers are unstable. The site- 
center dark breather is unstable for low enough frequen- 
cies and linearly stable with an even potential near the 
edge of the phonon band. 

At last, we have analyzed small amplitude bright and 
dark breathers superposed to uniformly stressed static 
states Xn = cn. In particular, the lowest frequency of 
bright breathers is inside the phonon band if < V"{c) < 
1, and the highest frequency of dark breathers is above 
the phonon band if V"(c) > 1. We have obtained local 
conditions on V for their existence as well as approximate 
analytical expressions. 
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